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Abstract 



The experimental ratio of ip' to J/ip is approximately a constant in pA 
collisions, but decreases as the transverse energy increases in nucleus-nucleus 
collisions. These peculiar features can be explained as arising from approxi- 
mately the same cc-baryon absorption cross section for ip' and J/ ip but greater 
disruption probabilities for ip' than for J/ip due to the interaction of the cc 
system with soft particles produced in baryon-baryon collisions. 



1 



Typeset using REVTgX 



High-energy heavy-ion collisions have become the focus of intense research because of 
the possibility of producing a deconfmed quark-gluon plasma during such collisions 
The J / if> suppression has been suggested as a way to probe the screening between a charm 
quark with its antiquark partner in the plasma ||. While the J /if) suppression has been 
observed [||,[| as predicted, the phenomenon can be explained by the absorption model ]6|-0 , 
which was actually introduced earlier |J to measure the total ip-N cross section using J /if) 
suppression. In the model of Gerschel and Hiifner ||, the collision of a J /if) particle with 
baryons of the colliding nuclei may lead to the breakup of the J /if) into an open-charm 
pair. An effective J/^-baryon absorption cross section cr^iifiN) of 5-6 mb can explain 
the systematics of the suppression phenomenon @,@]- One can alternatively describe J /if) 
suppression in terms of the interaction of the J / if) particle with produced hadrons (comovers) 
|7|.|8|]. A comparison of the production of if)' with J /if) has been suggested to distinguish 
between deconfmement and absorption ||10|| . 

The NA38 experimental measurements using protons and heavy ions at 200A GeV and 
450A GeV reveal three features JTTIJT^] : 1) if)' /if) is approximately a constant in pA collisions, 



independent of energies, 2) if)' /if) decreases as the transverse energy Ej> increases in SU 
collisions, and 3) if>'/if> for SU collisions is about 0.5 of that for pA collisions. The first feature, 
further supported by other pA experiments [ T3| , implies that in pA collisions if)' is suppressed 



in the same way as if). We would like to describe an absorption model with a^^ip'N) = 
c&bsiipN) and with additional soft-particle disruption to explain all three features of the 
phenomenon. 

The production of J /if) or if)' occurs by the interaction of the partons of one baryon with 
the partons of the other baryon. The incipient cc pair is created with a radial dimension 
of the order of ~ 0.06 fm at t c5 . It is necessary for the incipient cc system to evolve to the 
bound state rms radius of 0.24 fm for if) at and 0.47 fm for if)' at |T^JT5[] . Because 



J /if> is produced predominantly in the central rapidity region |16j, the incipient cc pair is 
produced predominantly in the central rapidity region. 

In soft particle production in a baryon-baryon collision, we envisage Bjorken's inside- 



outside cascade picture |L7j or Webber's picture of gluon branching |L8| as a q and a q 
(or a diquark) pull apart. After the collision, the diquark of one nucleon and the valence 
quark of the other nucleon pull apart and the gauge field between them is polarized. Gluons 
are emitted at t g , and the system hadronizes at th- The characteristics of these produced 
gluons cannot yet be determined from nonperturbative QCD, but the shape of the rapidity 
distribution of produced gluons should be close to that of the produced hadrons. Thus, 
produced gluons are found predominantly in the central rapidity region. Because we shall 
use the J/ip (or if}') production rate in a nucleon- nucleon collision as a unit of reference, it 
is not necessary to include explicitly the interaction of the incipient cc system with gluons 
and their hadronized products in the same nucleon-nucleon collision, when we study pA and 
nucleus-nucleus (AB) collisions. 

The space-time diagram for a typical pA collision is depicted schematically in Fig. la. 
The trajectory of an incipient cc pair, which is produced predominantly in the central 
rapidity region of the colliding baryons, does not cross the trajectories of soft particles 
produced in earlier or later collisions. Therefore, there is little interaction between the 
produced cc system and these soft particles. However, the cc system collides with baryons 
crossing its trajectory to lead to the breakup of the cc system into DDX. Such a reaction 
requires the production of at least one light-quark pair and is an inelastic process. In the 
collision at 200A GeV, the cc rapidities are separated from the baryon rapidities by about two 
units and the reaction cross section can be calculated in the Additive Quark Model (AQM) 
T9| , whose approximate validity and connection with soft Pomeron exchanges have been 



re-assessed recently from the systematics of total hadron- and 7-hadron cross sections |20 
Using the Glauber theory and a Gaussian thickness function, the total cc-baryon inelastic 
cross section in the AQM is given by Eq. (12.27) of Ref. 0: 



Cabs 



[cc- N) = -2* P 2 j:( Q )(-fT 7n, (1) 

n=l W 



where / = a cq / [3 2 = a cq /2'ir(f3% d + (3% + j3 2 q ) , a cq is the inelastic cross section for the collision 
of c (or c) and a constituent q of the baryon, a/3/3 C c and VSPn ar e the rms radii of the 



cc and the baryon respectively, and (3 cq is the c-q interaction range. We find below that 
c r abs(V'-^) = 4.2 mb. Taking V3(3n = 0.74 fm [^T| and neglecting (3 cqi we obtain from Eq. 



(0) ° cq = 0.753 mb, which leads to ^^(ip'N) = 4.27 mb, and cr a bs(cc — N) = 4.17 mb for the 
initial cc at V3Pcc = 0.06 fm. Thus, the absorption cross section is approximate the same 
for any cc state during all stages of its evolution because 6a cq << 2n(3 2 . Consequently, in 
pA collisions, if)' is suppressed in the same way as ip and ip'/ip is a constant independent of 
A and collision energies, in agreement with experimental observations. 

The approximate equality of the absorption cross section for different cc states is sup- 
ported by the experimental ratio (Jtotai{'4 } 'N)/atota\('4 , N) ~ 0.75 to 0.86±0.15, for ^/s ranging 



from 6.4 GeV to 21.7 GeV, as deduced from the photoproduction of vector mesons |22| , |23 
The approximate equality of the absorption cross section for the incipient cc and other cc 
states implies that a small incipient cc system is not transparent to the hadron medium, in 
agreement with the absence of color transparency for small hadron systems as indicated by 
experimental data in A(e,e'p) reactions at high Q 2 [p4j| . 



To study AB collisions, we adopt a row-on-row picture and consider a typical row with 
a cross section of the size of the nucleon-nucleon inelastic cross section, Oi n = 29.4 mb. 
The space-time diagram of the collision can be depicted schematically in Fig. lb. The 
trajectories of the cc system cross the trajectories of colliding baryons, and the process of 
absorption due to the cc-baryon interaction is the same in pA as in AB collisions. However, 
there are collisions, such as the ones at E and F in Fig. 16, where the trajectories of 
incipient cc systems produced there cross the trajectories of the produced soft particles. It 
is necessary to consider the additional interaction of the cc systems with soft particles in AB 
collisions but not in pA collisions. Because the cc systems and the gluons (or their hadronized 
products) are produced predominantly in the central rapidity region, their rapidities are not 
much separated and their relative kinetic energies are not large. At these low energies, the 
absorption of a produced gluon by the cc system, the screening of c from c by gluons, and 
the inelastic gluon scattering which excites the cc system to higher levels will contribute to 
the breakup of the cc system. As gluons carry color, the cross sections for these breakup 



processes increase with the color dipole moment of cc which is proportional to the separation 
between c and c. Furthermore, the threshold for if)' breakup is small compared to those for 
ip and Xi.2- Thus, the breakup probability for a cc(if>') system due to cc-gluon interactions 
at low energies is greater than those for J / if) and x- 

The hadronized product (comovers) of the produced gluons can interact with J/ip and 
if)' to lead to their breakup The breakup of if)', xi,2 and J /if) into DD require the 

threshold energy of 52, ~ 200, and 640 MeV respectively. In cc-hadron interactions at low 
energies below thresholds, J /if) and x cannot be broken up by low energy pions. Above the 
thresholds, the interaction is mediated by color gluon exchange which probes the color dipole 
moment of the cc system. Thus, the breakup probability due to cc-hadron interactions at 
low energies is also larger for the if)' system than those for the J/ if) and the x systems. This 
is different from the higher-energy (Pomeron-exchange dominated) Glauber case discussed 
earlier where if)' and J/ if) have about the same inelastic baryon cross sections. 

Assuming straight-line space-time trajectories and uniform distribution of baryons in 
nuclei, we find 



7^i = / a\ b m) { l ~ [ l - T ^ h >^ N )] A }{ 1 - [^-T B (b-b A )^ ba (m] B }F(bA), (2) 
where Tj^{bA) is the thickness function of A and the disruption factor F{b^) is 

i JV< n n 

F(b A ) = ——J2a(n)Y,eM-0 E + (3) 

iV > iV < n=l »=1 j=l,j^i 

Here, Nyi^b^) and N^b^) are the greater and the smaller of the (rounded-off) numbers of 
target nucleons ATa^a)*?™ and projectile nucleons BTs(b — b^)^ in the row at 6^ with 
the cross section of a in , and a(n) is obtained by simple counting to be 

a{n) = 2 for n = 1, 2, iV< - 1, and a(iV<) = A^> - iV< + 1 . (4) 

In Eq. (^|), we have assumed that when a cc(ip) system is produced in the j-th collision and 
soft particles are produced at the same spatial location in the z-th collision, the bound state 
survival probability is related to the cc-gluon interaction time tfj and the cc-hadron inter- 
action time tfj by ^^a 1 ^^^ 1 ^ where k^ g and are rate constants, averaged over the 



interaction history of the cc(if>) system. The interaction times, which must be nonnegative, 
are 



U + t h - Max(tj + tg, tj + t c5 ) , 



(5) 



t% = t n + t f - Max(*< + t h , tj + t c5 ) , 



(6) 



where ti = t± + (i — l)2m^d/ y/s, is the nucleon mass, d(= 1.93 fm) is the inter-nucleon 
separation in a nucleus, \fs the nucleon-nucleon center-of-mass energy, and tf is the freeze- 
out time. In Eq. (§), the step function 9 = Q(AT A (b A )a in - l)Q(BT B (b- b A )a m - 1)Q(A- 
1)Q(B — 1) is introduced to insure that there is no soft-particle disruption in pA collisions. 



if) into if)' • F° r simplicity, we do not treat xi,2 separately so that the extracted parameters 



It is not yet possible to ascertain the exact nature of all suppression mechanisms in 
AB collisions because of the uncertainties in the reaction cross sections (see below) and 
the characteristics of produced gluons. Besides the cc-baryon absorption, the suppression 
can be attributed to (A) produced gluons, (B) both produced gluons and hadrons, (C) 
produced hadrons (as in the comover model @||), or deconfined matter with no baryon 
absorption ||26|| . In our model, with (Jab s {ipN) = cr abs ('0 / A^) =4.2 mb fixed by pA data and 
a set of plausible time parameters t g = 0.1, th = 1.2, tf = 3, and t c5 = 0.06 (in units of 
fm/c), we obtain results calculated with rate constants (in c/fm) (A) k^ g = 0.2, k^ g = 3, 
k^ h = k^ h = (gluon disruption only), (B) k^ g = 0.2, k^ g = 1, k^ h = 0, k^ h = 1 
(gluon disruption for J /if) but gluon and hadron disruption for if>'), and (C) k^ g = k^ g = 0, 
k^h = 0.12, k^'h = 3 (hadron disruption only). For the case when all impact parameters 
are summed over, the quantity Ba^/AB is plotted as a function of A 1 / 3 + B 1 ^ 3 in Fig. 
2 for the three cases considered. The presence of additional soft-particle disruption in AB 
collisions relative to pA collisions is consistent with the experimental data in Fig. 2. To study 
if)'/ if) data for SU collisions at 200A GeV, we relate the transverse energy Ej, approximately 



The expressions for the production of if)' 



can be obtained from Eqs. (0-^) above by changing 



are actually for a "J/ if)" system with the observed if) : x = 62 : 30 mixture |25| . 
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to the impact parameter as given in Refs. |28] and |TT|, and the ratio ip'/ip is corrected 



for the feeding of J/tp by ip' . The theoretical results for B'a(ip')/Ba(ip) as a function of 
Ej. for the cases of (A), (B), and (C) differ by about one percent and are represented for 
simplicity by a single solid curve in Fig. 3. The theoretical ratio decreases with increasing 
Ej, and coincides with the pA limit, in good agreement with data. When we include all 
impact parameters, we obtain theoretically [B'a{ip') /Ba{J /ip)] su /[B'a(ip') /Ba(J /ip)] pA = 
0.62 for the three sets of parameters, approximately consistent with the experimental ratio 
[B'a(^)/Ba(J/^)] su /[B'a(f)/Ba( J/^)} pA = 0.52 ± 0.07 gl|. 



The parameter sets in (A)-(C) suggest greater disruption for ip' than J ftp due to their 
interaction with soft-particles. To resolve the ambiguities, it is interesting to note that while 
heavy-quark production by hadron-hadron collisions is inhibited by the OZI rule, there is no 
such inhibition in gluon-gluon collisions. The fusion of energetic gluons produced in different 
baryon-baryon collisions can lead to additional charm and strangeness production |2£| and 



may explain the enhanced charm and dilepton production in AB collisions relative to pA 
collisions observed in pt|. 

For a medium m whose constituents move with an average velocity v m with respect to 
cc(ipi), the rate constants k^ iTn is p m v m a^ im where p m = (dN m / dy) /7i\Rq£ is the contribution 
to the medium number density from a single NN collision and cr$ irn is the ipi-m breakup 
cross section. The time t Q is the mean point in the time when the medium exists in the 
specified form; t = (t g + i/J/2 = 0.65 fm/c for gluons, and t = (th + t/)/2 = 2.1 fm/c 
for hadrons. We can estimate p m by taking dN^/dy = 2dN g /dy = 2.3 and Rq = 0.5 fm. 



We can estimate v m by assuming that a produced gluon has a mass M ~ 1 GeV f29fl . At 



a temperature T = 200 MeV which is often found in these nuclear collisions, the most 



probable velocity is v m = ^2kT/M = 0.6 for gluons and v m ~ 1 for hadrons. Then, the rate 
constants k^. m suggest approximate order-of-magnitudes of a^ g ~ 1.4 mb, o^ g ~ 20 mb for 
case (A), a^ g ~ 1.4 mb, a^i g ~ a^h ~ 7 mb for case (B), and a^h ~ 0.9 mb, a^u ~ 21 mb 
for case (C). The excessively large if)' cross sections required to explain the if)' suppression 
in cases (A) and (C) may make the scenario (B) tentatively a more attractive description. 



The geometrical model |yj predicts the total hadron-A cross section at high energies 
proportional to the rms hadron radius and suggests (J^s^ip'N) much larger than cXabs^A), 
which is also assumed in the comover model §]. Using PQCD, Kharzeev and Satz pEfl 
claim that ipN total cross section for \fs = 6 GeV is about 0.3 mb. The geometrical model 
results and PQCD results differ from those obtained here. A recent calculation using an 
exchange potential gives a J/^-hadron dissociation cross section about 7 mb at 0.8 GeV 
kinetic energy |32| , in variance with the PQCD results of Ref . |26j . Much work remains to 
be done to resolve the differences. 



ACKNOWLEDGMENTS 

The author would like to thank T. C. Awes, F. E. Barnes, and Chun Wa Wong for helpful 
discussions and Dr. C. Lourengo for valuable comments. This research was supported by the 
Division of Nuclear Physics, U.S. D. O. E. under Contract DE-AC05-84OR21400 managed 
by Lockheed Martin Energy Systems. 



S 



REFERENCES 

[1] Quark Matter '95, edited by A. M. Poskanzer, J. W. Harris, and L. S. Schroeder, 
published in Nucl. Phys. A590 (1995). 

[2] For an introduction, see e.g. C. Y. Wong, Introduction to High-Energy Heavy-Ion Col- 
lisions, World Scientific Publishing Company, 1994. 

[3] T. Matsui and H. Satz, Phys. Lett. B178, 416 (1986). 

[4] C. Baglin et al, NA38 Collaboration, Phys. Lett. B220, 471 (1989); C. Baglin et al, 
NA38 Collaboration, Phys. Lett. B272, 449 (1991). 

[5] M. C. Abreu et al., NA38 Collaboration, Nucl. Phys. A544, 209c (1992). 

[6] C. Gerschel and J. Hiifner, Phys. Lett. B207, 253 (1988); C. Gerschel and J. Hiifner, 
Nucl. Phys. A544, 513c (1992). 

[7] R. Vogt, M. Prakash, P. Koch, and T. H. Hansson, Phys. Lett. B207, 263 (1988); S. 
Gavin and R. Vogt, Nucl. Phys. B345, 104 (1990); R. Vogt, S. J. Brodsky, and P. Hoyer, 
Nucl. Phys. B360, 67 (1991); R. Vogt, Nucl. Phys. A544, 615c (1992). 

[8] S. Gavin, Nucl. Phys. A566, 383c (1994); S. Gavin, H. Satz, R. L. Thews, and R. Vogt, 
Zeit. Phys. C61, 351 (1994). 

[9] R. L. Anderson et al, Phys. Rev. Lett. 38, 263 (1977); see also K. Gottfried and D. R. 
Yennie, Phys. Rev. 182, 1595 (1969). 

[10] S. Gupta and H. Satz, Phys. Lett. B283 439 (1992). 

[11] C. Lourenco, Proc. of the Hirschegg '95 Workshop, Hirschegg, Austria, 1995, CERN 
Report CERN-PPE/95-72, 1995 (LIP Preprint 95-03, 1995). 

[12] C. Baglin et al., NA38 Collaboration, Phys. Lett. B345 617 (1995). 

[13] D. M. Aide et al, E772 Collaboration, Phys. Rev. Lett. 66, 133 (1991); A. G. Clark 



9 



et al, Nucl. Phys. B142, 29 (1978); L. Antoniazzi et al, Phys. Rev. D46, 4828 (1976); 
H. D. Snyder et al, Phys. Rev. Lett. 36, 1415 (1976); K. J. Anderson et al, Phys. Rev. 
Lett. 42, 944 (1979); M. H. Schub et al, FERMILAB-PUB-95-058 (1995). 

[14] E. Eichten et al, Phys. Rev. D21, 203 (1980). 

[15] J.-P. Blaizot and J.-Y. Ollitrault, Phys. Lett. B199, 499 (1987). 

[16] N. S. Craigie, Phys. Rep. 47, 1 (1978), Table 1. 

[17] J. D. Bjorken, Lectures presented in the 1973 Proceedings of the Summer Institute on 
Particle Physics, edited by Zipt, SLAC-167 (1973). 

[18] B. R. Webber, Nucl. Phys. B238 492 (1984). 

[19] E. M. Levin and L. L. Frankfurt, JETP Lett. 2, 65 (1965); H. J. Lipkin and F. Scheck, 
Phy. Rev. Lett. 16, 71 (1966). 

[20] A. Donnachie and P. V. Landshoff, Phys. Lett. B296, 227 (1992); A. Donnachie and P. 
V. Landshoff, Phys. Lett. B348, 213 (1995); H. J. Lipkin, Phy. Lett. B335, 500 (1994). 

[21] T. T. Chou, Phys. Rev. D19, 3327 (1980). 

[22] M. Brinkley et al, Phys. Rev. Lett. 50, 302 (1983). 

[23] C. W. Wong and C. Y. Wong, (to be published). 

[24] T. G. O'Neil et al, Phys. Lett. B351, 87 (1995). 

[25] L. Antoniazzi et al, Phys. Rev. Lett. 70, 383 (1993). 

[26] D. Kharzeev and H. Satz, Phys. Lett. B334, 155 (1994); H. Satz, CERN Report CERN- 
TH/95-24 1995, |hep-ph 9502321 . 



[27] J. Badier et al, Zeit. Phys. C20, 101 (1983). 

[28] C. Baglin et al., NA38 Collaboration, Phys. Lett. B251, 472 (1990). 



10 



[29] C. Y. Wong and Zhong-Qi Wang, ORNL Preprint ORNL-CTP-95-06/hep-ph-9510233, 
1995. 

[30] C. Lourengo et al, NA38 Collaboration, Nucl. Phys. A566, 77c (1994); M. A. Mazzoni, 
HELIOS-3 Collaboration, Nucl. Phys. A566, 95c (1994). 

[31] B. Povh and J. Hiifner, Phys. Rev. Lett. 57, 1612 (1987). 

[32] K. Martin, D. Blaschke, and E. Quack, Phys. Rev. C51, 2723 (1995). 



11 



FIGURES 

FIG. 1. Schematic space-time diagram in the nucleoli- nucleon center-of-mass system, with the 
time axis pointing upward, (a) is for a pA collision and (6) is for an AB collision. The trajectories 
of the baryons are given as solid lines, the trajectories of an incipient cc system produced in some 
of the collisions are represented by thick dashed lines and the trajectories of soft particles produced 
in some of the baryon-baryon collisions by thin dashed lines. 

FIG. 2. The quantity Bajf^/AB as a function of A 1 / 3 + B 1 / 3 for pA and AB collisions. The 



data points are from the NA3 Collaboration |27| and the NA38 Collaboration ||. The solid curve 



gives theoretical results for pA collisions. For AB collisions, the theoretical results are shown as 
the long-dashed curve for cases (A) and (B) and as the dotted curve for case (C). 

FIG. 3. The ratio B' 'cr{ip') /Bcr(tfj) as a function of the transverse energy in SU collisions at 



200A GeV. Data points are from Ref. [11]. The theoretical results are shown as the solid curve. 
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